Neutrino Telescopes as a Direct Probe of Supersymmetry Breaking 
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We consider supersymmetric models where the scale of supersymmetry breaking lies between 5 
x 10 6 GeV and 5 x 10 s GeV. In this class of theories, which includes models of gauge mediated 
supersymmetry breaking, the lightest supersymmetric particle is the gravitino. The next to lightest 
supersymmetric particle is typically a long lived charged slepton with a lifetime between a microsec- 
ond and a second, depending on its mass. Collisions of high energy neutrinos with nucleons in the 
earth can result in the production of a pair of these sleptons. Their very high boost means they 
typically decay outside the earth. We investigate the production of these particles by the diffuse 
flux of high energy neutrinos, and the potential for their observation in large ice or water Cerenkov 
detectors. The relatively small cross-section for the production of supersymmetric particles is par- 
tially compensated for by the very long range of heavy particles. The signal in the detector consists 
of two parallel charged tracks emerging from the earth about 100 meters apart, with very little 
background. A detailed calculation using the Waxman-Bahcall limit on the neutrino flux and real- 
istic spectra shows that km 3 experiments could see as many as 4 events a year. We conclude that 
neutrino telescopes will complement collider searches in the determination of the supersymmetry 
breaking scale, and may even give the first evidence for supersymmetry at the weak scale. 
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Introduction — The origin of the radiative stability of 
the weak scale is one of the most important questions in 
particle physics today. A natural answer requires new 
physics at the TeV scale. Among the candidate theories, 
weak scale supersymmetry remains the most attractive 
scenario. Although this is in no small measure due to its 
theoretical appeal (it is a simple and natural extension 
of the usual space-time symmetries), it is also favored 
by data from electroweak observables. These point to 
a weakly coupled Higgs sector, one without significant 
deviations from the standard model in regard to elec- 
troweak precision observables. However, supersymmetry 
must be broken since the superpartners have not yet been 
observed. The supersymmetric spectrum is determined 
by the supersymmetry breaking mechanism. 

Supersymmetric models typically have a symmetry, 
called R-parity, which ensures that the Lightest Super- 
symmetric Particle, the 'LSP', is stable. Which of the 
supersymmetric particles is the LSP? This is determined 
by the scale of supersymmetry breaking, which we denote 
by \Z~F, and which can lie anywhere between 10 3 GeV 
and 10 12 GeV. If supersymmetry is broken at high scales 
such that y/F is larger than 10 10 GeV the LSP is typically 
the neutralino. If however supersymmetry is broken at 
lower scales, \/~F < 10 10 GeV, the LSP is typically the 
gravitino. In models where the LSP is the gravitino, the 
Next to Lightest Supersymmetric Particle (NLSP) tends 
to be a charged slepton, typically the right-handed stau. 
Since the NLSP decays to gravitinos through interactions 
that are suppressed by powers of V~F, if the supersymme- 
try breaking scale is high its lifetime can be quite large. 
In gauge-mediated SUSY breaking, for instance, we have 

/~F V /100 GeV\ 5 
mf R J 



10 7 GeV 



10 km , (1) 



where rrif R is the stau mass. Thus, for VF > 10 7 GeV 
if these NLSPs were to be produced by very high energy 
collisions they could travel very long distances before de- 
caying. In the last several years many interesting and 
realistic scenarios have been proposed in which the scale 
of supersymmetry breaking y/~F is low and could lie be- 
tween 5 x 10 6 GeV and about 5 x 10 8 GeV. These include 
models of gauge mediation^]], and gauge and Yukawa 
mediation[2|, warped higher dimensional models in which 
supersymmetry is broken on an infrared brane and there- 
fore the scale V~F has been warped down 0, theory space 
realizations of higher dimensional models y| and models 
of supersoft supersymmetry breaking which are charac- 
terized by Dirac gauginos 5j . In all these classes of models 
the NLSP is typically a right handed stau. 

The existence of diffuse fluxes of high energy neutri- 
nos, possibly associated with the production of cosmic 
rays, has been widely discussed in the literature. Col- 
lisions of these high energy neutrinos with nucleons in 
the earth at energies above threshold for supersymmetric 
production frequently result in the production of a pair 
of supersymmetric particles, which promptly decay into 
NLSPs. These typically have a high boost 7nlsp — 1000 
or larger and therefore will not decay inside the earth pro- 
vided the supersymmetry breaking scale y/F > 10 7 GeV. 
For 5 x 10 6 < \/F < 10 7 GeV a significant fraction of 
the decays will occur inside the earth. Since the NLSP 
is charged, its upward going tracks could in principle be 
detected in large ice or water Cerenkov detectors, such 
as ICECUBE 6]. This is in analogy with the standard 
model charged current interaction giving muons, the pri- 
mary signal in neutrino telescopes. 

Naively, one would expect that the event rates for 
NLSPs -which typically have masses in the 100 GeV 
range- are negligible when compared to those for muons, 
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since their production cross section must be considerably 
smaller than that of the standard model interactions. 
The reason is that the SM interactions come primarily 
from very small values of x, the parton momentum frac- 
tion, whereas the supersymmetric process is limited to 
x > m 2 /2MpE v , with to given essentially by the sum 
of the produced supersymmetric particles. This naive 
expectation is however, misleading. The crucial observa- 
tion is that the range of a slepton is much larger than 
that of a muon, since energy loss due to radiation sets 
in at much higher energies. In neutrino telescopes, muon 
events must be produced either right outside the detector 
or in it, since the muon range for the energies of interest 
is in the few to tens of kilometers. Thus, most of the up- 
going CC events produced in the earth are lost. On the 
other hand, the range of NLSPs is typically in the hun- 
dreds to thousands of kilometers. Then, unlike for the 
muon 0, a significant fraction of the NLSPs produced 
will range into the detector. 

In what follows we compute the number of NLSP 
events. For this purpose we calculate the supersymmetric 
production cross sections and the NLSP range. 

The SUSY Cross Section — In the scenarios under con- 
sideration, every vN interaction producing supersym- 
metric particles will result in a pair of NLSPs, which have 
a very long lifetime. In what follows, we will assume this 
lifetime to be large enough so that NLSPs do not decay 
in the earth. For simplicity, we also neglect mixing with 
Higgsinos in the gaugino sector. The dominant process is 
analogous to the SM CC interactions and corresponds to 
the t-channel exchange of charginos producing i>N — > i q, 
as shown in Figure ^a) and (b). The neutrino, always 
produced left-handed by the weak interactions, can in- 
teract either with a left-handed down- type quark (a), or 
with a right-handed up- type quark (b). This results in 
the partonic cross sections: 
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where s, t and u are the usual Mandelstam variables, 
and Myj, m| and to^ are the chargino, the left-handed 
slepton and the squark masses respectively. The left- 
handed slepton and the squark decay promptly to the 
lighter "right-handed" slepton plus non-supersymmetric 
particles. We also include the subdominant neutralino 
exchange, Figure H ( c )-(d). We take m«j = 250 GeV, 
mg = 250 GeV and three values for the squark masses: 
niq = 300, 600 and 900 GeV. These are very represen- 
tative values in the scenarios under consideration. Typ- 
ically, the tr is the NLSP, being heavier only than the 
ultra-light and very weakly coupled gravitino. Charginos 
and neutralinos tend to be heavier since they also feel 
the SU(2)l interactions. Finally, squarks are heavier 
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FIG. 1: Feynman diagrams for supersymmetric particle pro- 
duction in vN collisions Charged current (chargino) interac- 
tions: (a) Left-left interaction requiring the insertion of the 
gaugino mass in the t-channel line, (b) Left-right interaction. 
Neutral current: (c), (d). There are analogous diagrams for 
anti-neutrinos as well as for strange and charm initial quarks. 
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FIG. 2: uN cross sections vs. the energy of the inci- 
dent neutrino. The curves correspond to m; = 250 GeV, 
mm — 250GeV; and for squark masses m, = 300 GeV (solid) 
, 600 GeV (dashed) and 900 GeV (dot-dashed). The top curve 
corresponds to the SM charged current interactions. 



still since their masses affected by the strong interactions. 
In Figure |21 we plot the cross sections for supersymmet- 
ric production in i/N interactions as a function of the 
neutrino energy. Also plotted for comparison is the SM 
charged current cross section. As advertised earlier, the 
SUSY cross sections are still suppressed with respect to 
the SM, even when well above threshold. 

The NLSP Range — Once produced by the vN inter- 
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actions in the earth, the NLSP pair should range into the 
detector, just as the muons produced by CC events 0. 
Charged particles lose energy due to ionization processes 
as well as through radiation. The average energy loss is 
given by || 

-^ = a(/? 7 )+c(/3 7 ) f3j , (4) 

here a and c characterize the ionization and radiation 
losses respectively, and are slowly varying functions of 
the energy. The ionization loss can be approximated by 

MeV cm 2 

a(/3 7 )~0.08 — — — (17 + 21n/3 7 ), (5) 
gr 

nd is rather independent of the particle mass. On the 
other hand, assuming c(/?7) ~ const., the radiative en- 
ergy loss can be written as c f3j = (b m )E. Thus 
bf/rrin = bf R m.f R , and the radiative energy loss for the 
NLSPs scales inversely with the mass. This results in 
a much larger range for the NLSP as compared to the 
muon. Current bounds on m,f R are just above 100 GeV. 
As a reference value we take rrif R — 150 GeV. Therefore, 
NLSPs produced hundreds, even thousands of kilometers 
away are within range of the detector. This is to be con- 
trasted with the fact that muons must be produced at 
distances not larger than tens of kilometers from the de- 
tector in order to be observed. As we will see, this will 
somewhat compensate the suppression of the SUSY cross 
sections observed earlier. 

Signals in Neutrino Telescopes — In order to compute 
the event rates in neutrino telescopes, we need to know 
the incoming neutrino flux. The presence of cosmic neu- 
trinos is expected on the basis of the existence of high 
energy cosmic rays. Several estimates of the neutrino 
flux are available in the literature. In most cases, it is 
expected that km 3 neutrino telescopes will measure this 
flux. Here, in order to present projections for the number 
of observed SUSY events, we make use of the Waxman- 
Bahcall (WB) limit || as an estimate of the cosmic neu- 
trino flux. We consider an initial flux containing both 
and v e (in a 2 : 1 ratio). Since the initial interactions (see 
Figure ^) produce and these are nearly degenerate in 
flavor, the flavor of the initial neutrino does not affect 
our results. For the same reason, the possibility of large 
mixing in the neutrino flux is also innocuous here. In 
order to correctly take into account the propagation of 
neutrinos and the NLSP Ir through the earth, we make 
use of a model of the earth density profile as detailed in 
Reference jjjj- 

In Figure [3] we show the energy distribution for the 
NLSP pair events for three choices of squark masses: 
300 GeV, 600 GeV and 900 GeV. Also shown are the 
neutrino flux at earth in the WB limit, as well as the en- 
ergy distribution of upgoing /x's. We see that, even for the 
heavier squarks, it is possible to obtain observable event 
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FIG. 3: Energy distribution of tr pair events per km 2 , per 
year. From top to bottom: m,q — 300, 600 and 900 GeV. 
Here, mf R =l50 GeV and rriw — 250 GeV. Also shown are the 
neutrino flux at earth and the u flux through the detector. In 
all cases we make use of the WB limit for the neutrino flux. 



TABLE I: Number of events per km 2 per year for the WB 
and MPR fluxes. The first column refers to upgoing muons. 
The last three columns correspond to upgoing NLSP pair 
events, for three different choices of squark masses: 300 GeV, 
600 GeV and 900 GeV. The number of muon events are given 
for energies above threshold for production of a 250 GeV £l 
plus a 300 GeV squark, ic, 1.6 x 10 5 GeV. 
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rates. In Table [I] we show the event rates for £r pair pro- 
duction per year and per km 2 . The rates are given for the 
WB flux as well as for the Mannhcim-Prothcroe-Rachen 
(MPR) flux [J{j> both for optically thin sources. For com- 
parison, we also show the rates of upgoing muons. Thus, 
km 3 Cerenkov detectors such as ICECUBE, appear to 
be sensitive to most of the parameter space of interest in 
scenarios with a relatively long lived NLSP. 

Since the NLSPs are produced in pairs very far from 
the detector and with a very large boost, typical signal 
events consist of two tracks separated by SR ~ L 8, with 
L the distance to the production point (~ 100—1000 Km) 
and 9 ~ PsuW/Pboost ~ 10~ 3 — 10~ 4 . If we consider L 
to be of the order of the NLSP range, then in the lin- 
ear regime 5R ~ constant ~ 100m. As we have seen in 
the discussion following eqn. JSJ, for very high energies 
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TABLE II: Number of events for extended ICECUBE [jj per 
year assuming v flux is given by the WB limit. The II and 
squark masses and the number of muons are as in Table Q 
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the range grows logarithmically with energy, leading to 
somewhat smaller values of 8R ~ (20 — 40)m. The track 
separation is then mildly sensitive to the stau injection 
energy [1^, Then, most NLSP events would consist of 
two parallel but well separated tracks, and are therefore 
expected to be very distinctive and different from back- 
grounds. In addition, we see from Figure|31that the signal 
events have a harder spectrum than the \x background. 

Conclusions — We have shown for the first time that 
neutrino telescopes are potentially sensitive to the rela- 
tively long-lived charged NLSPs which are present in a 
wide variety of models of supersymmetry breaking. The 
event rates shown in Tableware already encouraging for 
experimental facilities that are been built, such as ICE- 
CUBE. The region of the supersymmetry breaking pa- 
rameter space that is available to neutrino telescopes is 
determined, on the one hand, by the twin requirements 
that the NLSP lifetime be long enough to give a signal 
(VF > 5 x 10 6 GeV), but not be so long as to disturb 
Big Bang Nucleosynthesis (VF < 5 x 10 8 GeV). Thus the 
observation of NLSP events at neutrino telescopes will 
constitute a direct probe of the scale of supersymmetry 
breaking. This is to be compared with the potential ob- 
servation of these NLSPs at the Large Hadron Collider 
(LHC), where for this range of y/~F, the NLSP decays 
outside the detector and is seen through its ionization 
tracks. However, the observation at the LHC would not 
constrain the NLSP lifetime significantly. Thus, we see 
that neutrino telescopes are complementary to collider 
searches. For instance, the observation of NLSP events 
at the LHC, coupled to the non-observation in neutrino 
telescopes would point to \/F < 10 7 GeV. 

Future upgrades of ICECUBE, as well as of the water 
detectors ANTARES Q and NESTOR [l^ will result 
in even better sensitivity. As an example, in Table ITT] we 
show rates for an expanded version of ICECUBE 13] . 

In the present letter we focused on supersymmetry. 
However, many other theories give rise to relatively long 
lived charg ed p articles which can be observed by neutrino 
telescopes [l2j. 

Acknowledgments — We would like to thank Azriel 
Goldschmidt and Ian Hinchliffe for useful conversations. 



I. A. is supported by NSF Grant Physics/Polar Programs 
0071886. G. B. and Z. C. are supported by the Direc- 
tor, Office of Science, Office of High Energy and Nu- 
clear Physics, of the U. S. Department of Energy under 
Contract DE-AC03-76SF00098 and by the N.S.F. under 
PHY-00-98840. 



[1] M. Dine, W. Fischler, M. Srednicki, Nuc. Phys. B 189, 
575 (1981); S. Dimopoulos, S. Raby, Nuc. Phys. B 192, 
353 (1981); L. Alvarez-Gaume, M. Claudson, M.B. Wise, 
Nuc. Phys. B 207, 96 (1982); M. Dine, A.E. Nelson, 
Phys. Rev. D 48, 1277 (1993); M. Dine, A.E. Nelson, 
Y. Shirman, Phys. Rev. D 51, 1362 (1995); M. Dine, 
A.E. Nelson, Y. Nir, Y. Shirman, Phys. Rev. D 53, 
2658 (1996). For a review, see G.F. Giudice, R. Rattazzi, 
Phys. Rept. 322, 419 (1999). 

[2] M. Dine, Y. Nir and Y. Shirman, Phys. Rev. D 55, 1501 

(1997) ; T. Han and R. J. Zhang, Phys. Lett. B 428, 
120 (1998); Z. Chacko and E. Ponton, Phys. Rev. D 66, 
095004 (2002). 

[3] T. Gherghetta and A. Pomarol, Nucl. Phys. B 586, 141 
(2000); W. D. Goldberger, Y. Nomura and D. R. Smith, 
Phys. Rev. D 67, 075021 (2003); Z. Chacko and 
E. Ponton, JHEP 0311, 024 (2003); Y. Nomura and 
D. R. Smith, Phys. Rev. D 68, 075003 (2003). 

[4] C. Csaki, J. Erlich, C. Grojean and G. D. Kribs, Phys. 
Rev. D 65, 015003 (2002); H. C. Cheng, D. E. Ka- 
plan, M. Schmaltz and W. Skiba, Phys. Lett. B 515, 
395 (2001); Z. Chacko, E. Katz and E. Perazzi, Phys. 
Rev. D 66, 095012 (2002); T. Kobayashi, N. Maru and 
K. Yoshioka, Eu r Phys. J. C 29, 277 (2003); N. Maru 
and K. Yoshioka, arXiv:hep-ph/0311337 

[5] P. J. Fox, A. E. Nelson and N. Weiner, JHEP 0208, 035 
(2002). 

[6] J. Ahrens et al. [The IceCube Collaboration], Nucl. Phys. 
Proc. Suppl. 118, 388 (2003). 

[7] For a detailed discussion of neutrino event rates see I. Al- 
buquerque, J. Lamoureux and G. F. Smoot, Astrophys. 
J. Suppl. 141, 195 (2002). 

[8] K. Hagiwara et al. [Particle Data Group Collaboration], 
Phys. Rev. D 66, 010001 (2002). 

[9] E. Waxman and J. N. Bahcall, Phys. Rev. D 59, 023002 
(1999); J. N. Bahcall and E. Waxman, Phys. Rev. D 64, 
023002 (2001). 

[10] K. Mannheim, R. J. Protheroe and J. P. Rachen, Phys. 
Rev. D 63, 023003 (2001) 

[11] R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, 
Astropart. Phys. 5, 81 (1996); R. Gandhi, C. Quigg, 
M. H. Reno and I. Sarcevic, Phys. Rev. D 58, 093009 

(1998) . 

[12] I. Albuquerque, G. Burdman and Z. Chacko, in prepara- 
tion. 

[13] F. Halzen and D. Hooper, arXiv:astro- ph/0310152 

[14] J. A. Aguilar et al. [T he~^NTARES Collaboration], 
arXiv:astro-ph/0310130 

[15] P. K. P. Grieder [NESTOR Collaboration], Nuovo Cim. 
24C (2001) 771. 



